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We present and investigate an extension of the classical random graph to a general class of 
inhomogeneous random graph models, where vertices come in different types, and the probability 
of realizing an edge depends on the types of its terminal vertices. This approach provides a general 
framework for the analysis of a large class of models. The generic phase structure is derived using 
generating function techniques, and relations to other classes of models are pointed out. 
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I. INTRODUCTION 

The concept of random graphs (RG) has recently be- 
come the target of an increasing interest, as a tool 
for modeling various kinds of networks, arising e.g. in 
physics, biology and biophysics as well as in social and 
information-technological structures. 

The classical RG model [1-3] describes a homogeneous, 
sparse random graph of order N, where each edge is ran- 
domly and independently realized with a fixed probabil- 
ity p — c/N. For large orders, there is a critical value of 
c = 1, above which almost every graph contains a single 
giant connected component being of order O(N), with 
the remaining components being small compared to N. 
This model yields an asymptotic degree distribution that 
is Poissonian with the average degree given by c. 

Many real-life networks, such as the internet, have been 
shown to possess other types of degree distribution, some- 
times displaying a power law behavior over many orders 
of magnitude, ruling out the classical RG as the relevant 
model. A number of alternative RG models have been 
suggested in an attempt to yield random graphs with 
more general types of degree distribution, such as the 
desired power behavior. Some of these models describe 
dynamical random graphs, where the graphs arise as the 
result of a stochastic growth process, such as randomly 
grown networks [4, 5], or scale- free networks based on 
preferential attachment [6]. Others focus on describing 
ensembles of random graphs with certain given proper- 
ties, without bothering about how they came about; a 
particularly interesting approach of this type, possessing 
a high degree of generality, is based on considering ran- 
dom graphs of fixed order with a given arbitrary degree 
distribution [7-10]. 

In this article, we will investigate a general class of 
sparse inhomogeneous RG models, by means of a straight- 
forward generalization of the classical model to a situa- 
tion where vertices may come in different types, such that 
the probability for an edge depends on the types of its 
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pair of terminal vertices. While this class of models inher- 
its certain features from the homogeneous model - such 
as the existence of a critical hypersurface in parameter 
space, beyond which asymptotically almost every graph 
has a giant component -, it is capable of producing a wide 
class of asymptotic degree distributions, among these dis- 
tributions with power law behavior. 

This general class of models is shown to contain a num- 
ber of existing models as special cases, and can be used 
as a general framework for the analysis of various RG 
models. 

The structure of this article is as follows. In Sec. 2, 
some of the more salient features of the classical model 
are briefly reviewed, while our generalization is presented 
and analysed in Sec. 3. In Sec. 4, a number of special 
cases are discussed, while Sec. 5 contains our conclusions. 



II. THE CLASSICAL MODEL 

Here will briefly review some of the more prominent 
properties of the classical random graph in the large N 
limit, to pave the ground for the subsequent analysis of 
its generalization. 

Definition 1 Let Q(N,c), with c a real positive number, 
denote the ensemble of graphs of order N > c, where each 
edge is independently realized with probability p = c/N . 

This ensemble has a critical value of c = 1, above which 
almost all graphs for large N have a single large con- 
nected component - the giant component - with a finite 
fraction of the vertices, while the remaining components 
are small. 



A. Exposing Connected Components 

The standard method to reveal the size distribution of 
the orders of components is to expose these as follows. 
Start with a single (random) vertex, reveal its neighbors 
by following edges, then their neighbors, etc. Let n& be 
the number of vertices exposed for the first time in step k 
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of this process. The distribution of rife, given the previous 
numbers n = 1, n 1; . . . rife_i, becomes 

P(n k )= ( N ~^l=O n l\ ( l _ q n k -^n k , n.-^N-Etom 

\ n k J 

(1) 

where q = 1 — c/7V. 

In the large TV limit with a fixed c, P(rifc) tends to 
e~ nk - lC (nk-ic) nk /rife!, and the process reduces to a Pois- 
sonian branching tree model <B(c), with each vertex inde- 
pendently branching to a number of new vertices, where 
this number is a Poissonian random variable with average 
c. The distribution p n over the order n of the resulting 
tree is conveniently analyzed in terms of the generating 
function F(z) = ^2 n p n z n , which must satisfy 

F(z) = zeMc(F(z)-l)). (2) 

This can be solved iteratively for each z, and F(z) must 
be a stable fixed point of the corresponding iterated map; 
it is easy to see that this implies < 1/ j cj . For 

\z\ < 1, there is a unique solution for F(z), reachable 
from the starting point 0, given by F(z) — C(z)/c, where 
C(z) is the unique solution to xe~ x = zce~ c in the unit 
disk. Expanding the corresponding inverse of xe~ x yields 
the exact result p n = n n ^ 1 c n ^ 1 e^ nc /n\. 

Particularly interesting is the result for z = 1, defin- 
ing / = F(l), which represents the total probability and 
might be expected to be 1, which is an obvious solution to 
Eq. (2) for z = 1. Indeed, for c < 1, this solution is sta- 
ble, but for c > 1, it becomes unstable, and another fixed 
point becomes the attractor, given by / = c/c, where c is 
the unique solution to ce~ c = ce~ c in the interval [0, 1]. 

Thus, for c < 1 the branching model is subcritical, and 
always terminates after a finite number of steps, while for 
c > 1 it is supercritical - the deficit in total probability 
is due to a finite probability 1 — / that the order n of the 
generated tree becomes infinite, i.e. that the branching 
process never terminates. 

For a large but finite N, this corresponds to all com- 
ponents being small, i.e. o(N), for c < 1, while for c > 1 
there exists a single giant component of order ~ N(l — f) 
with the remaining components being small, having an 
order distribution similar to that obtained for the com- 
plementary c-value c. 

I 



P{ni,k) = 




where qij = 1 — Pij = 1 — Cij/N w exp(— c^/iV). This 
expression can be simplified in different domains. 



III. GENERALIZATION TO 
INHOMOGENEOUS RANDOM GRAPHS 

The classical RG model can be generalized in a 
straight-forward way to inhomogeneous graphs by as- 
suming that vertices can come in different types i <G 
{1, . . . , K}. This enables us to consider a very general 
class of inhomogeneous RG models, to be referred to as 
IRG: 



Definition 2 Given a positive integer K , a K- 
dimensional vector r = {ri . . . rx} of positive probabil- 
ities summing to 1, and a symmetric K x K matrix c 
with non-negative elements Cij, let G(N, K, r, c) denote 
the ensemble of graphs G of order N, defined as follows: 

i) Each vertex is independently assigned a type i G 
{1 . . . K} with probability ri. 

ii) Independently for each unordered pair of vertices, the 
corresponding undirected edge is realized with probability 
Pij = c^/N, where (i,j) is the corresponding pair of ver- 
tex types. 



Remark 1 An asymptotically equivalent alternative is 
to fix the number of vertices of each type to certain 
values Ni « Nri, and possibly also the number of 
edges between vertices of types i,j to fixed values Eij w 
(l-SijWdjNiNj/N. 



A. Revelation of a Connected Component 

In analogy to the classical Q (N, c) model, the model 
Q(N , K,r,c) can be analyzed by recursively revealing a 
connected component by exploring neighbors, starting 
from a single vertex. Let n^fe be the number of new ver- 
tices of type i revealed in the fcth stage of the revelation 
(so rii.o = Si.i , with iq the type of the starting vertex). 
Given the number of revealed vertices of different types 
in the previous stages, n^k obeys the conditional distri- 
bution 



(3) 

I 

B. Small Components and the Branching Process 
Approximation 
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As long as the order of the revealed part is small as 
compared to N, we can approximate Eq. (3) by the Pois- 
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son distribution 

P(ni,k)«e-SiW-.iJ / . ( 4 ) 

This corresponds to approximating the revelation process 
by a Poissonian random branching process. 

For the distribution P(n) of the total number n of gen- 
erated vertices stating from a random vertex, we can de- 
fine a generating function F(z) = J2 n P(n)z n . Since the 
distribution will depend on the type of the initial vertex, 
F must be written as the weighted average of the corre- 
sponding generating functions Fi(z) for the distributions 
Pi{n) conditional on the initial type i, i.e. 



F(z) = J2r i F i (z), with Fi(z) =J2 P i( n > 



(5) 



in absolute value; the graphs then a.a.s. possess no giant 
component. 

B) When some eigenvalue is larger than one, the model is 
supercritical, and the graphs a.a.s. possess a giant com- 
ponent; its number Hi of vertices of type i asymptotically 
satisfies m/N ~ r,(l — fi), where the fi correspond to a 
stable solution of Eq. (7). 

Remark 3 It appears natural to require in addition that 
c cannot be block- diagonalized; otherwise ergodicity would 
be broken, and the graph would trivially decompose into 
distinct subgraphs, which could be treated separately. 

In the supercritical case, the generating functions Fi{z) 
can be renormalized with fi, to yield generating functions 
for the finite (non-giant) component part. Let Fi(z) = 
Fi(z)/fi. Then F is a stable solution of 



The vector F having the different F t as components sat- 
isfies the coupled set of equations 



Fi(z) = zexp X/vO/J/O - 1) 



(8) 



Ft(z) = zexp J^cyrj {F 3 {z) 



which is the inhomogeneous version of Eq. (2). 



(6) 



Remark 2 A more elaborate set of generating functions 
Fi{zx . . . zk) — J2n Pi( n ) Ylj ^j 3 could be defined, using 
a distinct variable z% for each type i, with rii the total 
number of revealed vertices of type i. These would obey 
equations obtained by replacing the "z" after the equal 
sign in Eq. ( 6) by '% ". Here we do not care about the 
detailed type content, and the simpler version, Fi(z) = 
Fi(z, . . . , z), will suffice. 

Interpreting Eq. (6) as a A'-dimensional iterated map 
(replace "=" by ":="), the proper solution is the sta- 
ble fixed point reached from the starting point F = 0. 
Particularly interesting is the result for z = 1, so let 
/, = Fi(l), expressing the probability that the branch- 
ing process will terminate, conditional on the type of the 
starting vertex, and let the unconditional counterpart be 
denoted by / = J2j r jfj = ^(1)- The fi satisfy the 
coupled set of equations 



fi = CX P CijTj {fj - 1) 



(7) 



with a naive solution f — 1, the stability of which can 
be analyzed by means of linearization of Eq. (7) around 
f = 1, yielding {c^rj} as the relevant matrix. This is all 
we need in order to pin down the appearance of the giant, 
as well as its asymptotic size, and we state the result 
without proof (it follows by analogy to the corresponding 
result for the classical model): 

Theorem 1 A) The model G(N, K, r, c) is subcritical if 
the eigenvalues of the matrix {cijfj} are all less than one 



with F(l) = 1. This describes a subcritical branching 
process with renormalized parameters fj = rj/,/r • f and 
Cij — c^ r • f . For a finite N, we must have N ~ N r • f , 
and we see that this conserves pij = Cij/N — p^; thus, 
the renormalized model is simply the naive restriction of 
the original one to the subset of vertices outside the giant 
component. 



C. Large Components and the Deterministic 
Approximation 

When the giant component is revealed, another ap- 
proximation can be made to Eq. (3). Once the number 
of revealed vertices become of O(N), the distribution of 
rii.fc becomes sharply peaked around its average, due to 
a factor of N appearing in the exponent. As a result, the 
fluctuations become negligible, yielding a deterministic 
iterative equation for the consecutive revealed numbers. 

In terms of the fraction gi ^ = 1 — X^=o n i,i/{^ r i) °f au 
vertices of type i not yet revealed after step k, this yields 



li,k = 9i,k-i cxp ^2 c ij r j(9j,k-i - 9j,k-2) , (9) 



revealing the conserved quantities 



9i,k exp c l3 rj (1 - g jt k-i) ■ (10) 



The values of /ii must be ~ 1, since their values can only 
change in an earlier stage when the number of revealed 
vertices is still small, but then g ~ 1; thus, in the large N 
limit we can safely assume fa = 1. The two-step recur- 
sion (9) reduces to a one-step recursion, taking the form 
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<7i,fc = e^j ■ Ci i r j^Sj,k-i i) ; wn i cn can be seen as iterating 
the map 

g % -» exp ^ CijTj- ( 5j - 1) j (11) 

until a stable fixed point is reached. If the model is sub- 
critical, this is given by the trivial fixed point gi = 1, 
whereas for a supercritical model a non-trivial fixed point 
with gi < 1 results, signalling the existence of a giant 
component containing a fraction 1 — gi of the vertices of 
type i. 

Eq. (11) is identical to Eq. (7), which was derived in 
the limit of small numbers of revealed vertices; thus, we 
have established the same set of equations in two different 
limits. 

Remark 4 A third, heuristic way of estimating the size 
of the giant component is as follows. Suppose the giant 
contains a fraction rii of the vertices of type i. Then we 
can estimate its neighborhood, i.e. the set of vertices con- 
nected to at least one vertex in the giant ( which of course 
must be the giant itself), as follows, based on the rather 
bold assumption that the edge probabilities do not depend 
on whether any or both of its terminal vertices are in the 
giant: The total number of vertices of type i is Nri . For 
each of these, the probability of not being connected to any 
of the vertices in the giant is cxp(— ^ . Cijnj/N). Thus 

we can expect a number Nn (l - exp(- J2 3 do n ] IN)) of 
vertices of type i in the neighborhood, i.e. in the giant. 
Writing n\ as Nn(l — gi), we recover Eq. (11), in spite 
of the bold assumptions involved. 

D. Extended Type Spaces 

While we have assumed a finite number of types K, 
defining the type space T = Zk, the above results should 
be more or less directly extendable to models where the 
type space T is a dcnumcrable infinite set, or even a 
continuous manifold, under some general conditions yet 
to be precisely determined. 

Definition 3 For a given type space T , with a normal- 
ized measure r on T , and a given non-negative symmetric 
function c(x , y) onT 2 , define Q{N ,T ,r,c) as the ensem- 
ble of RGs of order N, where each vertex is independently 
assigned a type x e T according to r(x), and for each ver- 
tex pair the corresponding edge is independently chosen 
with probability c(x,y)/N, with (x,y) the corresponding 
pair of types. 

For the denumerable case, T = Z + , it appears natural to 
require the asymptotic degree averages , or at least the 
total averages C, = YlJLi c ij r j > to be uniformly bounded. 
For cases where the elements of c are unbounded, an 
alternative is to regularize pij for finite N by using pij = 
1 — exp(— Cij/N) instead of the unbounded Cij/N. 



Also, reasonable care may have to be taken that c is 
sufficiently ergodic. Let tij be if = 0, 1 other- 
wise. The matrix t then describes a graph in type space, 
with tij = 1 corresponding to the existence of the edge 
(i,j). Then, sufficient ergodicity could e.g. mean that 
this graph should have a finite diameter, i.e. a uniformly 
bounded distance between vertex pairs. 

For the case of a continuous type space T, similar care 
must be taken. In addition, some kind of continuity con- 
straint seems appropriate, both on c and r. 

Note that a continuous T allows for a continuous repa- 
rameterization invariance. Thus, for the case of T = R, 
assume / to be a strictly increasing, continuously diffcr- 
entiable mapping of R to itself. Then the model defined 
by c(x,y) = c(f(x)J(y)) and r(x) = r(f(x))f'(x) is 
completely equivalent to the one with c(x,y) and r(x). 
Thus, r(x) could be transformed to any desirable normal- 
ized distribution on R. In particular, it could be trans- 
formed to the uniform distribution on the unit interval, 
yielding a kind of standard representation of the model. 
For higher-dimensional manifolds, things arc more com- 
plicated, and it appears difficult to devise a universal 
standardization procedure. 

A precise determination of feasibility conditions for ex- 
tended type spaces will be the subject of future work. 

E. Degree Distributions 

Many properties (but not all!) of a graph ensemble are 
reflected in its asymptotic degree distribution. In IRG, 
the asymptotic degree distribution p m is determined by 
r and c, and given simply as the weighted average of the 
type-specific degree distributions p m \i, being Poissonian 
with an average Cj defined by Cj = c ij r j- The result 
is 

Pm = Y / r i exp(-C i )^- (12) 
ml 

with the associated generating function 

H(z) = Y / PmZ m = Y^r i exp(C i (z-l)). (13) 

m i 

This puts a limitation on the possible degree distribu- 
tions that can be obtained within IRG: It must be pos- 
sible to write the distribution as a positive linear combi- 
nation of Poissonians, i.e. 

Pm = —. I c m e- c p(c)dc, (14) 
ml J 

where p(c) > describes an, a priori arbitrary, distribu- 
tion of type-specific Poissonian degree averages c = Cj, 
assuming the possibility of a continuum of types. This is 
a kind of smoothness constraint. In particular, it implies 

2 TO + 1 . . 

P m < Pm+lPm-l (15) 
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for each m > 0. While this excludes e.g. random regular 
graphs where the degree is fixed, it does allow for a wide 
class of degree distributions, such as distributions with a 
power tail, p m oc m~ a for large to, by letting p(c) having 
a similar power tail, p(c) oc cT a for large c. 

Note that a particular model in IRG is not determined 
solely by the degree sequence, which depends on Cij only 
through the average Cj = ^ • CijTj. This is in contrast to 
a class of recently considered models [7-10]; such models 
define a particular subclass of IRG, however, as will be 
shown below. 



IV. SPECIAL CASES OF INTEREST 

For K = 1 , of course the known properties of the clas- 
sical RG model is recovered. Below we will consider a 
few less trivial examples. 

A. Random Bipartite Graph 

Assuming two distinct vertex types, i.e. K = 2, a 
simple ensemble of random bipartite graphs results from 
the choice 



c = 



a 
a 



(16) 



With an arbitrary choice of type distribution r = (n, r 2 ), 
this yields for the asymptotic generating function F(z) = 
(Fi,F 2 ) the equations 



Fi(z) 

For z — 1, this yields 



zexp(ar 2 (F 2 (z)-l)), (17a) 
zexp(ari (Fi(z) - 1)) . (17b) 



fx = cxp (ar 2 (/ 2 - 1)) , h = cxp (an (/i - 1)) , 

(18) 

yielding the critical value of a as a c = l/^/n^- In the 
symmetric case of n = r 2 = 1/2, we have a c = 2, and 
/l = h = f satisfying / = cxp (a/2 (/ - 1)). 

In a similar way, ensembles of random if-partite 
graphs can be defined, which can be seen as generated by 
the complete graph JCk ■ Similarly, ensembles of random 
graphs based on an arbitrary generating graph can be 
defined, with c proportional to the incidence matrix for 
the generating graph. A nice twist results from using a 
random graph as a generator. 



Rank-1 c Matrix 



A particularly interesting special case results when c 
has the factorized form Cjj = CiCj/C, where d > can 
be interpreted as a connection tendency for vertices of 
type i, while C = J2i r iGi- 



Writing the asymptotic generating function as F(z) — 
^\ nFi(z), we get for Fi(z) in this case 



Fi(z) = zexp 



CiEj^CjjFjjz)-!) 

c 



(19) 



which can be reduced to a single equation for the function 

G(z) = ^ i r i C i F i (z)/C, reading 



G{z) = z J2 nCi exp (C t (G(z) - 1)) /C. 



(20) 



In terms of the generating function H(z) for the asymp- 
totic degree distribution, Eq. (13), this can be written 
as 



G(z) 



H (G(z)) 



(21) 



ff'(l) ' 
and in terms of G(z) we have 

F(z) = zH(G(z)). (22) 
For z = 1 in particular, we get for g = G(l) the equation 



.9 = 



H'( 9 ) 

H'(iy 



(23) 



and linearization around the trivial solution g = 1 yields 
stability for H (1)/H (1) < 1, corresponding to the 
model being subcritical for (C 2 ) < (C), which is equiva- 
lent to (to 2 ) < (2m) in terms of moments of the degree 
distribution. 

With c restricted to have rank 1, the resulting models 
are asymptotically equivalent to models from a superfi- 
cially very different class of random graphs that has re- 
cently attracted some attention [7-10]. There, a random 
graph ensemble based on an arbitrary asymptotic degree 
distribution p m is defined for a finite order N by ran- 
domly selecting a member from the set of graphs with a 
given degree sequence, such that the number of vertices 
with degree m is approximately Np m . 

Also for such a model, the recursive exposition of a con- 
nected component asymptotically yields a well-defined 
branching process, apparently very different from the 
Poissonian ones obtained for IRG. Here, the inital ver- 
tex is assigned a random degree to according to p m , and 
subsequently branches to m daughter vertices. Each new 
vertex is independently assigned a degree n > 0, dis- 
tributed according to np n / J2 m mp m (consistent with the 
assumption that the asymptotic probability of connect- 
ing to a particular vertex is proportional to its degree), 
and then branches to n — 1 daughters (since one of its 
edges is already used). 

The asymptotic generating function F(z) = J2k PkZ k 
for the resulting order distribution Pk then satisfies the 
equation 



F(z) = zJ2PmG(zY 



(24) 
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expressing the choice of the initial degree m. Here, G(z) 
is the edge generating function, which satisfies 

G(z) = z E "!? mG(z)m " (25) 

expressing the choice of the daughter's degree m, and its 
branching to m — 1 edges. 

These are nothing but Eqs. (22) and (21) in disguise, 
showing the complete asymptotic equivalence of the two 
models, despite the superficial differences; indeed, the 
criteria (m(m — 2)) < for subcriticality derived above 
arc in complete accordance with the results of Ref . [8] . 

C. Dynamical Random Graph with Finite Memory 

The last example is given by a recently proposed class 
of dynamical random graphs [5] with memory, where a 
graph is produced starting from a single node according 
to the combination of three random processes in contin- 
uous time, all Poissonian: 

1. For each existing vertex, new, initially isolated, ver- 
tices are added at a rate 7. 

2. For each existing vertex, new random edges are 
added at a rate A, connecting it to random existing 
vertices. 

3. Each existing edge is deleted at a rate [i. 

It is easy to see that the expected order of the graph 
grows with time t as e 7 *, and after an initial transient, 
vertices are only distinguished by their age, and we are 
led to consider a inhomogeneous model with a continuum 
of vertex types, T = [0, 00 [, given by vertex age. 

1. Asymptotic properties 

The probability density for ages x is asymptotically 
given by 

r(x) = 1 e-~< x . (26) 

For each pair of vertices, the probability of a connection 
is independent of the existence of other connections, and 
depends on the age x of the youngest vertex involved, and 
amounts to, at time t, p(x) = ( e (7-^)z _ 1 ) e -7*. 
We obtain 

c{x, y) = ( e^-ri min(x ' y) - l) , (27) 

7 - fj, \ J 

which seems feasible enough: c is ergodic, continuous and 
although c(x, y) is not uniformly bounded, the average 
C{x) = J c(x,y)r(y)dy = ^(1 - e"^) is (at least for 
yu > 0). Thus, we are lead to consider the spectrum of 
the integral kernel 

G(x, y) = f 1 - e-^-T) ™< x -y)\ e -™, (28) 

u — 1 \ > 



which is recognized as being proportional to the Green's 
function (i.e. a kernel representation of the formal op- 
erator inverse) for a particular differential operator C on 
R+, given by 

£ — 5^s(s + "-*)- (29 » 

with boundary conditions /(0) = 0, and /(a^e^ -7 ^/ 2 
growing at most as a power of x as x — > 00. Criticality 
results when the ground state of C has eigenvalue 1. 

With a finite memory, /i > 0, the eigenvalue equation 
for £ is a disguised version of Bessel's equation of order 
J //J,— 1 in the variable y = -\/8A7e~ iiX / 2 //i, and criticality 
results when the first positive zero X 7 / /J _ 1 of J 7 / M _i is 
given by V8A7/V, i- e - for A = M^^-i/ 8 !- 

In the special case of infinite memory, fi = 0, the model 
reduces to a randomly grown network [4], and yields 

£ =-*MK^H- <30) 

with eigenfunctions of the form e"* x / 2 sin{ojx) with eigen- 
value 7 , yielding the ground state value ^ for 
lu = 0, and criticality for A = 7/8. 

The above results arc all consistent with those obtained 
in Rcfs. [4, 5] on the phase structure for these models. 



V. CONCLUSION 

We have investigated a generalization of the classical 
homogeneous model of sparse random graphs, obtained 
by imposing a type structure on the vertices. This yields 
a very general class of inhomogeneous random graph 
models, and the asymptotic degree distributions are not 
restricted to Poissonians, but allow for various types of 
behaviour, within certain limitations. Thus, e.g., power 
behavior is possible, while a fixed degree (regular graph) 
is ruled out. 

The models in this class are not determined by the de- 
gree distribution alone, but contains an infinity of models 
for each possible distribution, in contrast to a recently 
considered class of models based on a given degree distri- 
bution. Interestingly enough, a relation does exist, since 
such models are shown to result in a special case of the 
present approach. 

In other special cases it describes the asymptotic static 
properties of certain models of evolving random graphs, 
such as randomly grown networks, and dynamical graphs 
with memory. 

Only certain aspects of the approach have been covered 
in this paper, and a more detailed analysis, e.g. of the 
feasibility conditions for extended type spaces, will be 
the subject of forthcoming work, as will investigations 
on further extensions of the approach e.g. to directed 
graphs. 
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